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Ultra-precise lattice parameter measurements in single crystals are achievable, in principle, by 
X-ray multiple diffraction (XMD) experiments. Tiny sample misalignments have hindered the sys- 
tematic usage of XMD in studies where accuracy is an important issue. In this work, theoretical 
basement and methods for correcting general misalignment errors are presented. As a practical 
demonstration, the induced strain of buried InAs quantum dots grown on GaAs (001) substrates is 
determined. Such a demonstration confirms the possibility to investigate epitaxial nanostructures 
via the strain field that they generate in the substrate crystalline lattice. 
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INTRODUCTION 



For decades, X-ray multiple diffraction (XMD) phe- 
nomenon has been applied in material science for study- 
ing a broad range of subjects: in-plane crystalline per- 
fection of surfaces and interfaces^; strain determination 
in single crystals caused by an external stimulus^ - —; the 
studies of habit modification of crystals due to the intro- 
duction of impurities species 7 ; direct observation of crys- 
tallographic phase transition in single crystals^; physical 
measurements of the invariant phase triplets^ - —; two- 
dimensional materials 14 or quasi-crystalsi£*i& and even 
in the plasma studies and diagnostics^ 7 ^. 

Among the reasons behind the success of the XMD in 
material science, one comes from the fact that XMD is 
the unique technique capable of providing information 
over several different directions within a crystalline lat- 
tice under the same diffraction geometry, i.e. with inci- 
dence ou and detector 26 fixed angles. Azimuthally scan- 
ning a chosen reflection, also known as Renninger scan- 
ning, allows accurate unit-cell lattice-parameters deter- 
mination. This information is very important for charac- 
terizing single crystals under extreme conditions, where 
the instrumental setup, crystal shape and/or limited 
choise of X-ray photon energies do not leave too many 
accessible reflections to be measured. For a clear ex- 
ample, in applying uniform-uniaxial electric fields, elec- 
tric contacts have to be set on two crystal faces. The 
field effects over the three-dimensional lattice can be 
studied by Renninger scanning a reflection that has not 
been hindered by the contacts^. Exploitation of the syn- 
chrotron radiation linear polarization to tune the rela- 
tive strength of the multiple diffracted waves has sig- 
nificantly increased the number of suitable reflections 
to be taken for Renninger scanning 19 . Moreover, it is 
also well known that this technique is the most accu- 
rate method for experimental lattice parameters mea- 
surements in single crystals22iSi; the precision can easily 



achieve 1/10000 A, depending on the chosen XMD case. 

High-precision lattice parameter measurements, abso- 
lute or relative (variational), are necessary in several re- 
search fields, such as determination of piezoelectric ten- 
sor, structural phase transition and dopping with low 
concentrations of impurities (typically parts per million) 
where the effects produced over the unit cell are very 
small. In studies where precise lattice parameter mea- 
surements are desired, it is essencial to account for all 
possible systematic errors in the Renninger scanning. A 
source of significant systematic errors originates from the 
misalingment between the azimuthal rotation axis 0-axis 
and the diffraction vector of the reflection to be scanned, 
named the primary reflection. Even in very good equip- 
ments and careful sample alignment procedures, residual 
tiny misalingments are always present. Since, in general, 
these misalignments can be smaller than 0.01°, they have 
been neglected so far and no further efforts were made 
to take them into account. However, as experimentally 
observed, even such small misalignment severely com- 
promise the accuracy of the results. In this article, we 
present the theory to account for a general sample mis- 
alignments as well as the experimental procedure for cor- 
recting systematic errors in lattice-parameter measure- 
ments by Renninger scanning. 



II. THEORY 

To excite multiple X-ray diffraction in crystals, the in- 
cident beam direction, wavevector fc, must fulfill two con- 
ditions summarized by 



and 



k P = -P P/2 



fe-S = -S-S/2. 
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P and S are the diffraction vectors of the primary and 
secondary reflections, respectively. Equations and 
J2J) stand for the Bragg cones of these reflections in the 
reciprocal space. In Renninger scanning, the primary 
reflection is kept excited, i.e. Eq. JU is fulfilled, while 
the crystal rotates around P; it means that the incident 
beam is parallel to the primary Bragg cone during the az- 
imuthal crystal rotation. The XMD occurs when the sec- 
ondary reflection is excited by the azimuthal rotation and 
therefore Eq. J2J) is also fulfilled. Note that more than 
one secondary reflection can be excited (3, 4,... beam 
cases) at the same azimuthal position; Eq. (J2J is valid for 
all simultaneously diffracting secondary reflections. 

The most well-known expressions for determining 
XMD positions in Renninger scans were obtained from 
Eq. (|2j)22i2£ assuming that the primary reflection is al- 
ways aligned, i.e. Eq. (QJ is fulfilled for a complete 360° 
rotation around P. However, it is not feasible to account 
for sample misalignments using the above couple of equa- 
tions. The reason is that Eq. is not sensitive to mis- 
alignments out of the incidence plane of the primary re- 
flection; in other words, when the misalignment between 
P and the the azimuthal rotation axis (</> axis) is taken 
into account by Eq. (QJ , it appears in Eq. Q as a negligi- 
ble second-order correction. Therefore, a third equation 
where the sensitivity to the three-dimensional misalign- 
ment of P is the same as in Eq. (|2| is necessary. Note 
that at least two equations are required since the XMD 
condition depends on the incidence uo and azimuthal <p 
angles i.e. depends on two variables. 

A third equation that is also fulfilled under XMD con- 
dition stands for the Bragg cone of the coupling re- 
flection, whose diffraction vector is given by C. Since 
P = S + C, it is possible to replace P by S + C in Eq. JQ) 
to obtain 




FIG. 1: Three-beam X-ray diffraction in misaligned crystals 
where the primary diffraction vector P is not perfectly aligned 
along the azimuthal scanning axis (0 axis). The rotation ma- 
trices Rx(Sx) and R y (S y ) account for the misalignment of the 
secondary and coupling diffraction vectors with respect to the 
(x, y, z) goniometer reference system, k is the wavevector of 
the incident radiation. 
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A common reference system is required to describe all 
vectors in the above equations. Here, the (x , y, z) go- 
niometer system shown in Fig. 1 will be used, z is along 
the <p axis, and x is an arbitrary reference direction cho- 
sen for <p = 0. For sake of simplicity, the vectors are 
represented by the 1x3 matrices of their components, 
for instance 



fc-C = -C-C/2-C-S. (3) 

In this equation, the Aco and A<fi deviations of the in- 
cident beam direction k{uj,(j)) have now about the same 
weight that they have in Eq. (0); it leads to a trivial sys- 
tem of linear equations for determining both deviations 
as a function of the sample misalignment. The depen- 
dence of the incident beam direction on these angular 
deviations can be taken into account by 



k ~ k + Aukuj + A0fe (4) 

where k u = dk/duo and k<j> = dk/d(j) are calculated 
at uoq and </>o, the incidence and azimuthal angles for 
exciting the XMD in non- misaligned samples, i.e. at 
fco = fc(cJo, </>o)- Any reciprocal vector G can also be 
written in terms of G°, its position in non- misaligned 
samples, and a rotation matrix R as G = G°R. Then, 
by using Eq. (0J into Eqs. (0) and JSJ, Auu and A(j) 
are determined by solving the following system of linear 



fc = — A 1 (cos uj cos 0, cos uj sin </>, sino;), (6a) 



k u = —A x (— sinc^o cos 0o, — sin ujq sin 0o ? cosc^o), 

(6b) 

fc^ = — A _1 (— cosc^o sin 0o 7 cos ujq cos 0o ? 0), (6c) 

and 

G° = (G x , G y , G z ). (6d) 

A general sample misalignment is taken into account 
by two rotation matrices R X (S X ) and R y (S y ) around the 
x and y directions (Fig. 1), and the total sample mis- 
alignment, for S Xi y <C 1, is 
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Therefore, 



AG = G-G° = (G^, G Z S X 



G X 5y 



GyS X ) 



(8) 



corresponds to the displacement of any diffraction vec- 
tor due to crystal misalignment. 

The incident beam direction for non-misaligned sam- 
ples fco = M^o^o) is calculated by Eqs. QJ and 
which provide uq and 0o since 



sinc^o = A|P| 2 /2 



and 



cos(^o — ol) = cos (3 = 



A|S| 2 /2-S^sincj 

S xy COS LdQ 



(9a) 



(9b) 



S xy cos a 



x • S° and S ?J 



sin a 



where 

v-s°. 

According to Eq. (|9b|) there are two azimuthal posi- 
tions, commonly called 'out-in' and 'in-out', where the 
same secondary reflection is excited: </>i = a — f3 and 
(j)2 = a + (3. Precise lattice-parameter determination pro- 
cedures are based on the experimental measurements of 
both positions. They provide 2/3 exp = <p2,exp — 4>i,exp for 
the calculation of the unit cell parameters from Eqs. (0, 
as described elsewhere (see for instance reference 20 ). The 
sample misalignment is measured by rocking curves of the 
primary reflection at (j) = 0, 90°, 180°, and 270°. If u , 
^90, ^180 and U270 are the respective rocking curves peak 
positions, 
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FIG. 2: (a) Renninger scans of the 002 GaAs reflection car- 
ried out with X-ray wavelength A = 1.330234A. The high- 
est peaks at 175.08° and 184.92° correspond to the excita- 
tion of the 111 ('in-out') and 111 ('out-in') secondary reflec- 
tions, respectively. The others secondary reflections shown 
in the diagram are also indexed. In the inset, [110] is the 
reference direction, = 0, and the vectors coresponding 
to the [111], [111], [111], and [111] directions are pointing 
outward the paper. The eight closed circles symbols in- 
dice the azimuthal (j) positions of all secondary reflections 
and their values are shown in Table 1. (b) Several high- 
resolutions cj) scans measured around <j) — 175.08° position 
(solid circle) and after a 360° (j) rotation (open circles). The 
peak positions have been determined by fitting the experi- 
mental intensity curves (open/solid circles) with lorentzian- 
gaussian convolution curves (solid lines), as explained in the 
text. The maximum shift observed in the peak position is 

§max - §min = 0.00075°. 



8x = (^270 - ^9o)/2 and S y = (uiso - wq)/2. (10) 

The Acj) n corrections in the 'out-in' (n = 1) and 'in- 
out' (n — 2) azimuthal positions are then obtained from 
Eq. (0), and they can be used either to refine the 2(3 exp 
value according to 

2Pexp = (4>2,exp ~ A0 2 ) - (<f>l,exp ~ A0i) (11) 

or to estimate misalignment effects on different sec- 
ondary reflections, for instance those providing some use- 
ful information about the crystalline structure of the 
sample. To demonstrate this in practice, the in-plane 
strain of a GaAs (001) substrate induced by In As quan- 
tum dots (QDs) is inspected in the following sections of 
this article. 



III. EXPERIMENTAL 

Data collection has been carried out at the Brazilian 
Synchrotron Light Laboratory (LNLS) D12A (XRD-1) 
beam-line, with the polarimeter-like diffractomei^. The 
smallest step size is 0.0004° in both uj and axes. The 



wavelength value, A = 1.330234 A from a Si (111) double- 
crystal monochromator, was measured by rocking curves 
the 111 and 333 Silicon reflections. The incidence plane 
was set to the vertical position, scattering upwards (cr- 
polarization). The sample used in this study is a comer- 
cial GaAs (001) substrate with InAs quantum dots (QD) 
grown by molecular beam epitaxy (MBE) using a rate of 
0.009 monolayer per second (mL/s). A 300A thick GaAs 
cap layer was grown on top of the QDs. 



IV. RESULTS AND DISCUSSIONS 

Table |D shows the 'out-in' (0i) and 'in-out' (02 ) az- 
imuthal XMD peaks positions for the 111, 111, 111, and 
111 secondary reflections in the Renninger scanning of 
the 002 GaAs reflection. The peak positions were de- 
termined by fitting the intensity data with lorentzian- 
gaussian convolution curves: a lorentzian standing for 
the intensity profile (FWHM = 0.0048°) of the peaks 
while the instrumental broadenig was taken into account 
by a gaussian function (FWHM = 0.0060°). A short por- 
tion of the Renninger scan, around = 180°, is shown 
in Fig. 2a; the [110] direction was taken as reference for 
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cf) = 0. To assure that the data are free of instrumental 
errors, the <ft scan of a chosen XMD peak, in this case the 
peak at <\> = 175.08°, was repeated after measuring the 
other peaks without changing the <ft axis rotation sense, 
i.e. after a 360° (j> rotation, see Fig.(2)b. Each <ft scan 
has been performed at the (002) rocking-curve maximum 
which was carried out about 0.5° before the XMD peak 
positions. 

The rocking curves used to characterize the sample 
misalignment are shown in Fig. 3; they provide S x ~ 
0.0009° and S y c± 0.0067°. However, in general, for small 
sample misalignments Eq. (0) can be linearized by nu- 
merical derivation, and the corrections of the XMD 
positions can be written as A(j) n = A n S x + B n 5 y . For the 
particular set of secondary reflections investigated here, 
the values of A n and B n coefficients are given in Table |I] 
(last two rows). By replacing these values into Eq. ([TT]). 
it is possible to demonstrate that 

Ano] = (Pin +Aii)/2 (12a) 

and 

Alio] = (0111 +Au)/2 (12b) 

do not depend on S x and S y . In other words, the 
average (3 exp for the secondary reflections with in-plane 
components in the [110] and [110] orthogonal directions, 
P[no] and /?[iio], are misalignment-free experimental val- 
ues, and therefore, very useful for precise lattice param- 
eter measurements in both orthogonal directions. 

We assume a tetragonal lattice distortion where the 
substrate in-plane lattice parameter is given by a = 
b = (1 — v)ao while, to preserve the unit cell vol- 
ume, c = (1 + 2v)a Q . For a = 5.6534A, |P| = 2/c, 
S° = (ft/a, fc/fo, 1/c) (in the crystal reciprocal space) 
and ft, k = ±1, Eq. Q provides 

$(y) = 85.08984° - 33.85V (13) 

Note that |S| = |S°| is invariant under any rotation, but 
their components are not. Therefore, S^, S y and S z in 
Eq. (|9b|) must be calculated for x = (1, 1, 0)/V% y = 
(1, 1, 0)/a/2 and z = (0, 0, 1), i.e., S x = (ft + fc)/(a>/2), 
S y = (-ft + k)/(aV2) and S z = 1/c. 

In the 002 Renninger scan of (001) samples, hkl sec- 
ondary reflections give rise to special XMD cases, named 
Bragg-surface diffractions, where the secondary-beam di- 
rections are very close to surface in-plane directions. The 
extreme asymmetry of such secondary reflections limits 
the secondary wavefield penetration depth to less than 
1 jim. Experimentally, this value has been measured 
by reciprocal-space mapping of the 002 rod in GaAs 
(001) crystals undergoing Bragg-surface diffraction. It 
has been demonstrated that about 80% of the diffraction 
intensity, for X-ray photons of 8 keV, is attenuated in the 
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FIG. 3: Experimental (open circles) and simulated (solid 
lines) rocking curves, o;-scans, of the 002 primary reflec- 
tion measured at <fi = 0, 90°, 180° and 270°, as shown 
aside of each scan. These scans were simulated by using 
two gaussian curves (dashed lines) whose splitting is about 
Acj = -0.00422(24)°. It provides Ac/c = 3.0(2) x 10" 4 . 
The misalignment angles calculated from the gaussians at 
lower u) values are S x ~ 0.0009° and S y ~ 0.0069°, while 
S' x (= 0.0024°) and S' y (= 0.0103°) were obtained from the p 
values in Table [I] 

first 0.3 iim bellow the surface 26 . According to this fact, 
the observed tetragonal deformation would correspond to 
an average value of a shallow layer just below the surface 
where the induced strain due to the InAs QD is signif- 
icant. Although the InAs epitaxial growth generates an 
expansive stress in the substrate lattice under the QDs, 
the adjacent regions are compressed as schematically il- 
lustrated in Fig. 21 

Differences in the tetragonal distortion along the [110] 
and [110] directions, i.e. /3[ 110 ] 7^ /%io]? would yield a 
unit cell sligthly twisted near the substrate surface. For 
instance, tiny variations in the 7 angle, between the a and 
b lattice vectors, can be estimated by numerical deriva- 
tion of Eq. (J2J. In this case, it leads to 

Alio] - Ano] = (C110 - Cii )A 7 (14) 

where Cx 10 = 0.1276 and Cno = —0.1302, and then to 
7 = 89.9985° ±0.0018°. 
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TABLE I: Azimuthal positions of the 111, 111, 111, and 111 secondary reflections in the Renninger scanning of the 002 GaAs 
primary reflection. Each position was measured three times (rows 1, 2 and 3) as explained in the text, cj) — (4>max + </>rmn)/2, 
£ = (4>max — <j)min)/ c 2<, ftexp = (02 — 0i)/2, fi are the average misalignment- free values, v stands for the unit-cell tetragonal 
distortion (Eq. C3J) and f3 m is = + (A0 2 - A0i)/2 (Eq. JTTJ), where A0 n = A n 5 f x + B n 6' y , S' x = 0.0024°, and S' y = 0.0103°. 
A n and B n were estimated by numerical derivation of Eq. Angular values are given in degrees. 
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FIG. 4: Schematic illustration of the compressive stress gen- 
erated by InAs QDs on GaAs substrate. There is also a 
300A thick GaAs cap layer over the QDs (not shown in the 
figure) . 



Finally, the agreement between the S x and S y and S' x 
and S f y values is better when rocking curves are measured 
with higher definition in their maximum peak position. 



V. CONCLUSIONS 



pie misalignment is eliminated and any systematic er- 
ror is averaged out. Therefore, a precision better than 
1/10000A can be easily achieved. 

As the method presented here can determine very small 
single crystal lattice deformations, it is an excellent tool 
for studying epitaxial nanostructures grown on top of 
crystalline substrates. The strain field produced by the 
nanostructures in the substrate lattice can be determined 
even when those nanostructures are buried by cap layers. 
As an application, QDs grown on GaAs (001) were an- 
alyzed and the results showed that the GaAs lattice is 
tetragonally distorted and slightly twisted near the sub- 
strate surface. 

The method however, is not restricted to the partic- 
ular geometry of the QDs system studied in this work 
where, due to the primary diffraction vector high sym- 
metry, a family of eight secondary reflections could be 
measured. When the system symmetry under study is 
lower, it is only important to obtain rocking- curves pro- 
files with well-defined peak positions. 



In summary, a method that corrects any sample mis- 
alignment and allows the determination of highly ac- 
curate single-crystals lattice parameters or their varia- 
tion was presented. The advantage of this method is in 
the measurement of rocking curves in orthogonal posi- 
tions. By combining those rocking curves and measure- 
ments of a XMD secondary reflection family, having the 
same (3 value and measured in their "in-out" and "out- 
in" azimuthal positions, the effects of a residual sam- 
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